Search for: All records

Abstract The Alpine‐Himalayan belt is one of Earth's most dynamic and complex regions, characterized by intense tectonic deformation and seismicity. Comprehensive analyses of continental‐scale crustal deformation and seismic hazards along this extensive orogenic belt require the compilation of large geodetic data sets. In this study, we integrate 42 published Global Navigation Satellite System (GNSS) velocity fields, building an internally consistent data set for the entire belt, spanning from Iberia to Southeast Asia and comprising 11,177 horizontal and 3,940 vertical velocities. We use this unified GNSS velocity field to estimate surface strain rates and their posterior uncertainties in the eastern Mediterranean region and the India‐Asia collision zone. Our results show large‐scale agreement between the orientation and style of geodetic and seismic strain rate tensors across the belt. Additionally, our analyses substantiate previously documented azimuthal alignments between principal strain rate directions and seismic anisotropy orientations, often used as a proxy for finite strain in the convecting mantle. These correlations are particularly apparent in the Aegean, North Anatolia, Tibet, Tian Shan, Altai, Sayan, and Baikal regions, underscoring the need for future research on the relationship between mantle flow and lithospheric deformation. 

Shear wave splitting of teleseismic core phases such as SKS is commonly used to constrain mantle seismic anisotropy, a proxy for convective deformation. In plate boundaries, sharp lateral variations of splitting measurements near transform faults are often linked to deformation within a lithospheric shear zone below, but potential seismic waveform effects from heterogeneous structure on small scales may influence the interpretation. Here, we explore possible finite frequency effects on shear wave splitting near fault zones in a fully three‐dimensional anisotropic setting. We find that shear zones wider than 80 km, a scale set by the Fresnel zone, can be clearly detected, but narrower zones are less distinguishable. Near the edge of the shear zone, the combined effect of anisotropy and scattering generates false splitting measurements with large delay times and fast axis orientation approaching the back‐azimuth, a bias which can only be identified when records from different back‐azimuths are analyzed together. This substantiates that back‐azimuthal variations of splitting can arise not just from vertical layering but also lateral changes of anisotropic media. We also test the effects of shear zone edge geometry, epicentral distance, filtering frequency, crustal thickness, and sediment cover. Our study delineates the ability of shear wave splitting to resolve and investigate fault zones, and emphasizes the importance of good azimuthal coverage to correctly interpret observed anisotropy. Based on revisiting previous shear wave splitting and lithospheric deformation studies, we infer that many crustal fault zones are underlain by lithospheric shear zones at least 20 km wide. 

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR. 

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper. 

Azimuthal variations in receiver function conversions can image lithospheric structural contrasts and anisotropic fabrics that together compose tectonic grain. We apply this method to data from EarthScope Transportable Array in Alaska and additional stations across the northern Cordillera. The best-resolved quantities are the strike and depth of dipping fabric contrasts or interfaces. We find a strong geographic gradient in such anomalies, with large amplitudes extending inboard from the present-day subduction margin, the Aleutian arc, and an influence of flat-slab subduction of the Yakutat microplate north of the Denali fault. An east–west band across interior Alaska shows low- amplitude crustal anomalies. Anomaly amplitudes correlate with structural intensity (density of aligned geological elements), but are the highest in areas of strong Cenozoic deformation, raising the question of an influence of current stress state. Imaged subsurface strikes show alignment with surface structures. We see concentric strikes around arc volcanoes implying dipping magmatic structures and fabric into the middle crust. Regions with present-day weaker deformation show lower anomaly amplitudes but structurally aligned strikes, suggesting pre-Cenozoic fabrics may have been overprinted or otherwise modified. We observe general coherence of the signal across the brittle-plastic transition. Imaged crustal fabrics are aligned with major faults and shear zones, whereas intrafault blocks show imaged strikes both parallel to and at high angles to major block-bounding faults. High-angle strikes are subparallel to neotectonic deformation, seismicity, fault lineaments, and prominent metallogenic belts, possibly due to overprinting and/or co-evolution with fault-parallel fabrics. We suggest that the underlying tectonic grain in the northern Cordillera is broadly distributed rather than strongly localized. Receiver functions thus reveal key information about the nature and continuity of tectonic fabrics at depth and can provide unique insights into the deformation history and distribution of regional strain in complex orogenic belts. 

Abstract The spatial distribution of the geochemical domains hosting recycled crust and primordial (high‐³He/⁴He) reservoirs, and how they are linked to mantle convection, are poorly understood. Two continent‐sized seismic anomalies located near the core‐mantle boundary—called the Large Low Shear Wave Velocity Provinces (LLSVPs)—are potential geochemical reservoir hosts. It has been suggested that high‐³He/⁴He hotspots are spatially confined to the LLSVPs, hotspots sampling recycled continental crust are associated with only one of the LLSVPs, and recycled continental crust shows no relationship with latitude. We reevaluate the links between LLSVPs and isotopic signatures of hotspot lavas using improved mantle flow models including plume conduit advection. While most hotspots with the highest‐³He/⁴He can indeed be traced to the LLSVP interiors, at least one high‐³He/⁴He hotspot, Yellowstone, is located outside of the LLSVPs. This suggests high‐³He/⁴He is not geographically confined to the LLSVPs. Instead, a positive correlation between hotspot buoyancy flux and maximum hotspot³He/⁴He suggests that it is plume dynamics (i.e., buoyancy), not geography, which determines whether a dense, deep, and possibly widespread high‐³He/⁴He reservoir is entrained. We also show that plume‐fed EM hotspots (enriched mantle, with low‐¹⁴³Nd/¹⁴⁴Nd), signaling recycled continental crust, are spatially linked to both LLSVPs, and located primarily in the southern hemisphere. Lastly, we confirm that hotspots sampling HIMU (“high‐μ,” or high²³⁸U/²⁰⁴Pb) domains are not spatially limited to the LLSVPs. These findings clarify and advance our understanding of deep mantle reservoir distributions, and we discuss how continental and oceanic crust subduction is consistent with the spatial decoupling of EM and HIMU.